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Рассматривается метод генерации одноуровневых марковских полей. Свойства порождаемого поля, выбор типа его составных 
элементов (алфавит) задаётся оператором. Предлагаемый алгоритм позволяет получать горизонтальные/вертикальные и диа- 
гональные типы мозаик. 


Введение 

При цифровом моделировании систем, пред- 
назначенных для обработки многомерных случай- 
ных полей (яркостные поля плоских изображений, 
температурные и другие поля атмосферы, океана и 
т. п. )ш , возникает необходимость машинной ими- 
тации тестовых сигналов в виде числовых полей с 
заданными спектральными и вероятностными 
свойствами. Сложность этих процессов такова, что 
их исследования аналитическими методами, по- 
зволяющими проникнуть в физику явлений, очень 
трудоёмки и не всегда могут быть доведены до кон- 
ца. Экспериментальные исследования слишком 
дороги, ограничены в возможности проведения 


факторного анализа (повторение эксперимента в 
строго определённых условиях просто невозмож- 
но) и в конечном итоге основаны на субъективных 
оценках ряда специалистов. Путь имитационного 
моделирования является одним из приемлемых 
объективных путей решения поставленной задачи 
в определённой предметной области. Он позволяет 
проводить быструю проверку гипотез, упрощать 
выкладки, а при создании макетов и планировании 
экспериментов позволяет проводить детальный 
анализ полученных результатов. Проведение фак- 
торного анализа, повтора экспериментов, при ис- 
пользовании цифрового моделирования не вызы- 
вает больших затруднений. 
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Рис. 1. Элементы упрощенного алфавита 



Известно, что реальным оптическим изображе- 
ниям соответствуют поля яркости, близкие к ку- 
сочно-постоянным функциям двух переменных, 
т.е. состоящие из областей, внутри которых яркость 
почти не меняется, и которые отделены друг от 
друга резкими границами-контурами [2]. При мо- 
делировании случайных полей применяют различ- 
ные способы их имитации. Разбиение производит- 
ся прямыми, геометрическими фигурами, с ис- 
пользованием фрактального разбиения (алгорит- 
мы Серпинского), посредством объединения опре- 
делённой общности точек на основе предиката од- 
нородности и т. д. Полученные области заполняют 
случайными постоянными яркостями. К основ- 
ным недостаткам этих подходов можно отнести 
сложность управления свойствами получаемых по- 
лей. В частности, сложно выделить векторную со- 
ставляющую контурного рисунка. Одним из воз- 
можных подходов к решению задачи получения те- 
стовых изображений является использование алго- 
ритмов, основанных на построении управляемого 
однородного марковского поля. Свойства порож- 
даемого поля связаны между собой формулами 
Пальма [3, 4], а выбор типа его составных элемен- 
тов (алфавита), а следовательно и типа решётки, 
производится оператором. Другими словами, ос- 
новной особенностью подобного типа алгоритмов 
является возможность управления вероятностны- 
ми и спектральными свойствами генерируемого 
поля (морфологией) [5, 6], получения векторного 
описания контурного рисунка, лежащего в его ос- 
нове. Необходимо отметить, что прежние алгорит- 
мы [5] не позволяли производить построение кон- 
турного рисунка с диагональными элементами. 
Этот факт затруднял использование алгоритмов ге- 
нерации марковских полей в проведении фактор- 
ного анализа определённых предметных областей. 
Например, при анализе детекторов контурных ри- 
сунков применение алгоритмов, описанных в ра- 
боте [5], не позволяет выявить закономерности де- 
тектирования диагональных элементов. 

Расширение алфавита, являющейся составной 
частью построения контурного рисунка, позволило 
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произвести создание нового по своим возможно- 
стям алгоритма, который успешно решает задачу 
формирования марковского поля не только с гори- 
зонтальными/вертикальными но и с диагональны- 
ми составляющими. 

Описание алгоритма 

Построение одноуровнего однородного марков- 
ского поля проводится с помощью ряда элементар- 
ных контурных элементов, называемых алфавитом 
построения. Полный алфавит состоит из 248 воз- 
можных элементов. Элементы выстраиваются в слу- 
чайную последовательность на основе системы 
условных вероятностей. Изменение вероятностей 
появления того или иного элемента алфавита позво- 
ляет генерировать различные типы контурных реше- 
ток от простейшей прямоугольной решетки до весь- 
ма замысловатых орнаментов. Для простоты демон- 
страции работы алгоритма рассмотрим упрощённый 
алфавит, состоящий из 23 элементов, рис. 1. 

Контурные элементы могут быть образованы 
лучами в восьми направлениях ( а , Ъ, с, сі, с,/, & к ), 
исходящих от центра 2 (рис. 2). 
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Рис. 2. Описание окрестности контурного элемента 

Величины а, Ь, с, сі, е, /, & к принимают значе- 
ния «О» или «1» и означают отсутствие или наличие 
соответствующего контурного луча. Для упроще- 
ния описание литер алфавита, удобно присвоить 
каждому литеру десятичный код. Значение деся- 
тичного кода определяется по формуле: 

N = а2* +Ъ2 1 +с2 2 +<12 ъ + 

+е2 4 + /2 5 + §2 6 + к2 7 . 

Свойства контурного рисунка (мозаики) пол- 
ностью определяются значениями финальных ве- 
роятностей при различных значениях а, Ь, с, сі, е , /, 
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& к. Для упрощения записи финальных вероятно- 
стей, бинарные значения а , Ъ , с, <і, е,/, & // приведе- 
ны ниже в десятичном виде. 

ДО) = г 0 ; 

/>(34) = />(136)=|; 

/>(10) = />(40) = />(130) = Д160) = 

/>(42) = />(138) = />(162) = />(168) = 

Р(170)=г 4 ; 

/>(21) = />(69) = />(81) = />(84) = 

/>(5) = />(20) = />(65) = />(80) = 

/>(17) = />(68) = у; 

/>(85) = г 8 ; 

8 

ІУ = Г 

і =0 

где г 0 - вероятность появления контурного элемен- 
та 0, г, - 34, 136, г 2 - 10, 40, 130, 160, г 3 - 42, 138, 162, 
168, г 4 - 170, г 5 - 21, 69, 81, 84, г 6 - 5, 20, 65, 80, 
г 7 - 17, 68, г 8 - 85. 

При построении контуров используется услов- 
ные вероятности Р(/§к/аЪсс1е ) , полученные из фор- 
мулы полной вероятности [7] 

Р(/ 8 к/аЪЫе ) = , 

XXX Р(аЬссІе/ 8 к ) 

/= 0 ^= 0^=0 


Полученные значения Р(/$к/аЬсс1е) являются 
элементами стохастической матрицы условных ве- 
роятностей Р размерностью 32x8. 

Построение контурного рисунка начинается с 
создания опорных исходящих литер в первой стро- 
ке и в первом столбце изображения. Левый верх- 
ний угол изображения считается началом отсчёта. 
Значение вектора (акссіе) для этого элемента не 
определено, а следовательно выбор литеры при его 
построении производится произвольно с учётом 
финальных вероятностей. На рис. 3 выделены се- 
рым цветом опорные задающие элементы, начало 
отсчёта обозначено крестом, стрелками обозначена 
координатная система. Генерация первой строки и 
первого столбца проводится по методике [5] с учё- 
том изотропности появления любого возможного 
контурного элемента. Вероятности появления кон- 
турных лучей р ю р ь , р п р й , р е или вероятности их от- 
сутствия д а , д ь < 7 С , вычисляются по формулам: 

Ча =Чс=Че=Л • 0 +г 1 +г 2 +г 3 +г 4 + ^-+Лі+^-; 

г. 2 к г, 

Чь=Чл =г 0 + 2 + ^ + 4Г +г 5 + г 6 + г 7 + г 8 ; 

Ра=' [ -Ча’Рь=' [ -Чь’ 

Рс= 1 ~Чс’Р 1 1 = 1 -Чл,Ре= 1 -Че- 

После создания первой опорной строки и пер- 
вого опорного столбца, производится построчное 
(слева на право, сверху вниз) прохождение узлов 
сопряжения контурных элементов. Исходя из за- 
дающего, априорно известного бинарного вектора 
(аЬЫе), определяется результирующий вектор (/§к). 
Эта процедура производится на основе условий, 
описанных в стохастической матрице условных ве- 







Рис. 3. Пример создания последовательности опорных элементов 
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роятностей. Заключающим этапом является зане- 
сение бинарной величины «2», значение которой 
определено полным бинарным вектором ( аЪЫе - 
/%к), в растр генерируемой мозаики. Рис. 4 поясня- 
ет геометрию проводимой генерации. 
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Рис. 4. Прямоугольный фрагмент изображения 


Последующие итерации повторяются по выше- 
определённой схеме с учётом сдвига прямоуголь- 
ного фрагмента в окрестности исследуемой обла- 
сти на одну точку в необходимом направлении. 

Достоинством предлагаемого алгоритма явля- 
ется возможность моделирования мозаичного изо- 
бражения с расширенным типом контурного ри- 
сунка. Эта возможность была достигнута за счёт 
ввода дополнительных элементов алфавита 


(рис. 1). Расширение базового алфавита позволило 
получать диагональные составляющие контурного 
поля, что является существенным достоинством 
предлагаемого алгоритма по сравнению с алгорит- 
мами, приведёнными в работах [5, 7]. Работоспо- 
собность алгоритма была проверена путем модели- 
рования мозаик, полученных в работах [5-7]. Мор- 
фология контуров задавалась исключением одного 
или нескольких типов контурных элементов в со- 
ответствии с методикой предложенной в работе [5]. 

Тип А является простейшей прямоугольной ре- 
шеткой, построенной на основе элементов {0, 170, 
34, 136}, рис. 1. Тип С образуется непересекающи- 
мися прямоугольниками разных размеров на осно- 
ве элементов {0, 34, 136, 42, 138, 162, 168}. Тип Д 
представляет собой многосвязные замкнутые обла- 
сти с более мелкими вкраплениями на основе эл- 
ементов {0, 34, 136, 0, 40 130, 160}. Мозаика /полу- 
чена на основе элементов {0, 17, 68}. 

Результаты моделирования 

Для моделирования различных типов мозаик 
был создан комплекс программ генерации одноу- 
ровневого марковского поля, написанный на алго- 
ритмическом языке С++ с использованием про- 
граммной среды моделирования Сапіаіа 2001. На 
рис. 5 приведены результаты моделирования кон- 
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Рис. 5. Разновидности мозаик 
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турных полей типа А, С, Л и К Контурные рисунки 
построены размерностью 200x200. Изображения 
сохранялись в формате *.Ътр. Заполнение незави- 
симыми яркостями произведено не было. 

Как видно из рисунков предложенный алго- 
ритм позволяет обобщить опыт работ [5-7], мозаи- 
ки А , С, Д а также получать однородное марков- 
ское поле не только с горизонтальными/верти- 
кальными но и с диагональными составляющими 
контурного рисунка, мозаика К 
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